Towards the fundamentals of car following theory 
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The problem of a car following a lead car driven with constant velocity is considered. To derive the 
governing equations for the following car dynamics a cost functional is constructed. This functional 
ranks the outcomes of different driving strategies, which applies to fairly general properties of 
the driver behavior. Assuming rational driver behavior, the existence of the Nash equilibrium 
is proved. Rational driving is defined by supposing that a driver corrects continuously the car 
motion to follow the optimal path minimizing the cost functional. The corresponding car-following 
dynamics is described quite generally by a boundary value problem based on the obtained extremal 
equations. Linearization of these equations around the stationary state results in a generalization 
of the widely used optimal velocity model. Under certain conditions (the "dense traffic" limit) the 
rational car dynamics comprises two stages, fast and slow. During the fast stage a driver eliminates 
the velocity difference between the cars, the subsequent slow stage optimizes the headway. In the 
"dense traffic" limit an effective Hamiltonian description is constructed. This allows a more detailed 
nonlinear analysis. Finally, the differences between rational and bounded rational driver behavior 
are discussed. The latter, in particular, justifies some basic assumptions used recently by the authors 
to construct a car-following model lying beyond the frameworks of rationality. 

PACS numbers: 89.40.-a, 45.70.Vn, 02.50.Le 

Keywords: car following model, motivated behavior, cost (utility) function, rationality, optimal driving 
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I. CAR-FOLLOWING THEORIES AND BASIC 
PROPERTIES OF DRIVER BEHAVIOR 

Recently, the theoretical and empirical foundations of 
the physics of traffic flow (for a review see Refs. 0, 0) 
has come into the focus of the physical community. The 
motion of individual cars has many peculiarities, since 
it is controlled by motivated driver behavior, together 
with some physical boundaries. Nevertheless, on macro- 
scopic scales the vehicle ensembles displays phenomena 
like phase formation and phase transitions widely met in 
physical systems (see, e.g., Refs. 0,0, 13)- So, the coop- 
erative behavior of cars treated as active particles seems 
to be of a more general nature than the mechanical laws 
and constructing a consistent theory of traffic flow "from 
scratch" is up to now a challenging problem. 

To describe individual car dynamics a great variety of 
microscopic models have been proposed. These models 
differ in the details of the interaction between cars and 
the time update rule, ranging from differential equations 
to cellular automata 0, U . There has been a big deal of 
work on the macroscopic behavior emerging from the mi- 
croscopic dynamics when exploring the behavior of sys- 
tems of interacting cars. However, there is still a lot of 
controversy in both the macroscopic behavior when com- 
pared to reality 0, and in the microscopic foundations 
of the individual car dynamics itself [5|- 

One currently adopted approach to specify the micro- 
scopic governing equations of the individual car motion 
is the so-called social force model. More details can be 
found in Refs. [f| 0, @; here only the basic ideas are 
touched on. At each moment t of time, a driver a changes 



the speed v a of her car depending on the road conditions 
and the arrangement of the neighboring cars: 
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The term g a {v a ) describes the motion of car a on the 
empty road, whereas the term g aa > (x a , v a \ x a i , v a > ) al- 
lows for the interaction of car a with car a' {a 1 ^ a). 
The interaction is due to the necessity for driver a to 
keep a certain safe headway distance between the cars. 
All the models mentioned above use various Ansdtze for 
the last term. 

The most interesting special case, which covers the 
majority of all traffic flow situations, is that of single- 
lane traffic. Here, all cars can be ordered according to 
their position on the road in the car motion direction 
x a < x Q +i, here a — 1, . . . , N. Most models take into 
account solely nearest neighboring cars a and a + 1, i.e., 
9aa' 7^ for a' = a + 1 and, may be, a' = a — 1 only. 
However, more complicated models exist that can be de- 
scribed as models with anticipation l^llfllllL Il2l IT^ | or 
the so-called intelligent driver model ElM" 

The earliest "follow-the-leader" modelsjla [T^ | relate 
the acceleration a a of car a to the velocity difference 
(v a - v a+ i) only, i.e., 



dv a 
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(1.2) 



where t v is the characteristic time scale of the veloc- 
ity relaxation. In subsequent generalizations of this 
model t v became a function of the car motion state, in 
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particular, of the current velocity v a and the headway 
h a = x a+ i —x a —i (for a review see Refs. 0,^1). Here, 
I is the car length. In Refs. [T^. l2fl | another approach 
called the optimal velocity model is proposed that de- 
scribes the individual car motion by 

~17~ = [Va~ ^opt(h a )} , (1.3) 

at t v 

where $ pt(M is the steady-state velocity (the optimal 
velocity) chosen by drivers as function of the headway h 
between the cars. It should be noted that this app roach is 
related to much earlier safety distance models |2lLl22ll23| . 

Concerning the fundamentals of approximations such 
as a a = a{v ai v a +i,x ai x a+ i), it is noted that there are 
actually two stimuli affecting the driver behavior. One 
of them is the necessity to move at the mean speed of 
traffic flow, in the given case at the speed v a +\ of car 
a + 1. So, first, driver a should control the velocity 
difference (v a — w Q +i). The other is the necessity to 
maintain a safe headway distance h opt (v a ) depending on 
the current velocity v a . The "following- the- leader" mod- 
els mainly take into account the former stimulus. The 
optimal velocity model, conversely, allows for the latter 
stimulus only . More so phis ticated approximations, e.g., 
Refs. [llllllllllllllllll^ to name but a few, allow 
for both stimuli. Note also a simple Ansatz called the 
combined model in Ref. J29] which is also related to the 
intelligent driver model |!4lll5|: 

dv a (1 - ft) r . 

—rr = [v a ~ v a+1 \ 

at t v 

- — [V a -1?opt(M] ' (!- 4a ) 

This equation takes into account both stimuli via a phc- 
nomenological coefficient < ft < 1. This is also the 
case for the Helly model [24| which can be written as (cf. 
Ref. i) 

dv a 1 

—rr = [ y a - v a+l\ 

at t v 

H t~ [(xa+i - x a ) - h op t(v a )] , (1.4b) 

T V L H 

where Lh is a certain spatial scale and h op t(v) is the 
optimal headway distance chosen by drivers when moving 
at speed v. Ref. |24| used a linear Ansatz for h opt (v). 
Later 26], this model was generalized to allow for the 
dependence of the kinetic coefficients on the motion state. 

However, the question whether a collection of variables 
such as {v ai Wa+i, x a , specifies the acceleration a a 

completely is not trivial. Drivers are characterized by 
the motivated behavior rather than physical regularities. 
For example, memory effects may be essential and can de- 
stroy the direct relationship a a — a(v a ,v a +i,x a ,x a +i). 
Up to now, memory effects in car following modelling 
have been treated only in a simplified version. This has 
been done by relating the current acceleration a a (t) to 



the velocities v a (t — r a ), v a+ i(t — r a ) and the headway 
distance h a (t — T a ) taken at a previous moment (t — r a ) of 
time (for a review of such an approach concerning with 
the "following-the-leader" models see Refs. [a fl8l| . r e- 
garding the optimal velocity model see Refs. |30l I3fll32fl ) . 
Here r a is the formal delay time in the driver response 
which is treated as a constant. Such an approach, how- 
ever, is rather formal, because, first, it is not clear why 
the memory effects relate only two moments of time in- 
stead of a certain interval as a whole. Second, a sim- 
ple physiological delay in the driver response as well as 
the mental driver estimation of the surrounding situa- 
tion should contribute to the value of r a . If the latter 
contribution is essential the delay time r a is to depend 
substantially on the state of car motion. Moreover, from 
our point of view the memory effects stem from the fact 
that the description of individual car motion is a bound- 
ary value problem rather than an initial value problem. 
Indeed, at the current moment of time the headway dis- 
tance and the car velocity are quantities given for the 
driver beforehand. Correcting the car motion she should 
choose such a driving strategy that in a certain time in- 
terval the car velocity and headway distance attain their 
optimal stationary values, at least, approximately. 

These memory effects are the topic of this paper. We 
propose that drivers plan their behavior for a certain time 
in advance (3j| instead of simply reacting to the sur- 
rounding situation. Similar ideas related to the optimum 
design of a distance controlling driver assistance system 
are discussed in Ref. [34]]. Mathematically, the driver's 
planning of the further motion is just to find extremals 
of a certain priority functional that ranks outcomes of 
different driving strategies. 

The derivation of microscopic governing equations for 
systems with motivated behavior based on a certain "op- 
timal self -org anizatio n" p rinciple has been discussed re- 
cently 0, Hi |1E EH!! H3- The assumption adopted 
in these works is that individuals try to minimize the in- 
teraction strength or, equivalently, to optimize their own 
success and to minimize the efforts required for this. The 
approach discussed here applies to the concepts of math- 
ematical economics, namel y, to the notion of p references 
and utility (see, e.g., Ref. |4£j). In Ref. 33] a specific 
form of the priority functional has been proposed. From 
that, a certain Ansatz like the combined model (|I.4al) or 
the Helly model l|1.4bjl could be derived. Nevertheless, 
the question how to find the priority functional "from 
scratch" remains open. 

In the present paper the priority functional is con- 
structed by applying to general properties of the driver 
behavior for a simplified situation. A car following (no 
overtaking allowed) a lead car which is driven with con- 
stant speed V is considered. The task is to derive gov- 
erning equations for the following car motion specified by 
the time dependence of the velocity v(t) and the headway 
distance h(t). 
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FIG. 1: Two equivalent car motion paths a and (3 illustrating 
the property of the cost functional with respect to path parti- 
tioning. Letters O and D stand for the origin and destination 
points, respectively, and the symbols {hi,Vi} denote the local 
states of traffic flow. 

II. THE COST FUNCTION OF DRIVING 

A. General properties of driver preference 

Assuming a driver to be able of comparing any two 
states {hi,vi}, {h 2 ,v 2 }, the phase plane {h > 0, v > 0} 
can be ordered by a preference relation ^. Therefore, the 
existence of a cost function Tih, v) may be assumed such 
that: 

{hi, vx} r< {h 2 , v 2 } T{h u v x ) > T{h 2 , v 2 ) . (2.1) 

Typically, ~-!F(h, v) is called the utility function and (for 
it) the relation -< matches the inequality <• Here, the 
use of condition (|2.1|) is preferred, because then the cost 
function T(h, v) is similar to the free energy of physical 
systems, with the minima as the stationary states. 

Obviously, the cost function T(h, v) cannot be speci- 
fied completely because a composite function ^ [J-(h, v)] 
also meets condition H2.lt f° r an Y increasing function 
W [■] . Therefore, at the current stage of the theory de- 
velopment any approximation or Ansatz adopted for the 
cost function T[h, v) has no meaning. 

To overcome this problem and to fix a certain class of 
the cost functions the preference relation defined on a set 
of car motion paths {h(t), v(t)} is considered. A car mov- 
ing from one point to another can meet different states 
of traffic flow. It enables one to unite various points on 
the phase plane in an element treated as a whole with re- 
spect to the driver preference, reducing the numerical un- 
certainty in their evaluation. These parts are supposed, 
first, to connect the same origin-destination pair {0,D} 
and, second, to go through a traffic flow pattern whose 
structure is known to the driver beforehand. Such a set 
obviously can be ordered by the driver preference and the 
corresponding cost functional £{h(t), v(t)} is assumed to 
exist. Here, only the fact that a car motion path can con- 
sist of various fragments of the phase plane is concerned 
with. Analysis of the real car dynamics is postponed to 
the following section. 

Figure^shows two paths a and f3 which are equivalent 
with respect to the driver preference. They differ from 



each other only in the order of passing through the traffic 
states, {hi, v±} — > {hi, Vi} — > {h 2 , v 2 } for the path a and 
{hi, vi} — ► {h 2 , v 2 } — ► {hi,Vi} for the path (3. It demon- 
strates the fact that if a path goes through a certain small 
neighborhood of the point {h, v} one or many times then 
only the cumulative time interval St during which the sys- 
tem was located in the given neighborhood contributes to 
the cost functional. The corresponding weight coefficient 
T{h, v) obviously plays the role of a certain cost function 
for the given traffic state. The resulting cost functional 
for the whole path is 

destination 

C{h(t),v{t)}= J «(*)]*• (2.2) 

origin 

If the driver has incomplete information about the traf- 
fic flow pattern ahead she plans the further motion only 
within some limits, spatial and temporal ones. This can 
be taken into account by introducing a certain cofactor 
in front of the function T(h, v) which depends on spatial 
coordinates or time. 

The initial introduction of the cost function remains 
the freedom of its choice too "wide" , as was noted above. 
The constructed form of the cost functional reduces this 
freedom substantially. Namely, let T(Ji, v) and !F{h, v) be 
two functions entering a cost functional similar to 12.2(1 
with the prediction effect such that they preserve the 
preference relation with respect to the car motion paths 
on an one-lane road. Then, as shown in Appendix ^ 
these cost functions are related by the expressions 

T{h,v) = Kj r (h,v)+C a (h)v, (2.3a) 
T{h,v) = KT{h,v)+C b {t), (2.3b) 

depending on the spatial or temporal limitation in the 
driver prediction. Here k > is an arbitrary posi- 
tive constant and C a (h), Cb(t) are some functions of the 
headway distance h or the time t. In other words, rela- 
tions l|2.3afl . i|2.3b|) specify the family of the admissible 
type of cost functions preserving the driver preference 
relation and, in addition, correspond to linear transfor- 
mation of the cost functional ( Appendix lA")l . 

For example, [£{h(t)}} 2 is also a cost functional. How- 
ever, formula (|2.2I) enables us to fix its specific form being 
linear with respect to time integration along car motion 
paths. We note that this construction conforms with re- 
garding the motion cost as the "physiological" cost of 
driving per time, i.e., having the meaning of a rate. 

Both of these families possess the same freedom in 
choosing a specific form of the cost function especially 
dealing with the extremals of the cost functional. In fact, 
the type (|2.3b|) is chosen in the following. 

B. Characteristic features of the car motion cost 

In the following, the priority function of the car mo- 
tion state with respect to the velocity v and the reciprocal 
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value p = 1/h of the headway distance will be character- 
ized. Note, that the value p is not the real car density on 
a highway. The present analysis ignores the car length, 
so we have preferred to use the value p defined as above. 
The most preferable state is the motion on an empty road 
(p = 0) at maximum speed $ m ax- This maximum speed is 
determined by external factors. So, at {p = 0, v = $ max } 
the cost function T (p,v) := J 7 {h,v)\ h _ x , has its global 
minimum: 

^(0,tf max ) =0, (2.4) 

set equal to zero keeping in mind the aforesaid about the 
freedom in specifying the cost function. Since there is no 
other minimum for the motion on an empty road, 

J"(0,0) = 1 (2.5) 

can be fixed. 

Because of the car construction a driver can visually 
control the headway distance within some value I ~ 1- 
2 m. When the headway distance h attains such values 
a driver has to stop her car because of possible collision 
even at sufficiently slow velocities. So, I coincides with 
the characteristic headway distance in dense jams where 
the car density attains the possible maximum. In the 
following, all headway distances are related to the scale 
I and the dimensionless variable pi is used. 

When not moving at all (v = 0) and the density of cars 
surrounding the given car does not come close to the limit 
values, i.e., pi <C 1, it does not matter how many cars are 
located in the vicinity. So the assumption that at v = 
and for h I the cost function is independent of p can 
be used: 

F{p, 0) = 1 for pi < 1 . (2.6) 

In all other cases, T (p, v) depends on certain combina- 
tions of p and v rather than on p individually, at least 
when pi <C 1. 

Considering the behavior of T (p, v) for a fixed speed 
v it can be stated that driving with small values of h 
(large values of p) requires a lot of effort. Thus, T (p, v) 
increases with p. The opposite case of small values of 
p (large values of h) deserves special attention. With- 
out the possibility of overtaking no especially attractive 
headway distance can be marked. Therefore, we assume 
that the cost function T (p, v) possesses the only one min- 
imum attained at the boundary point ,0 = provided the 
velocity v is fixed. 

Keeping in mind the aforesaid it is reasonable to write 

T(p,v)=T(0,v) + f(-^)(pir (2.7) 

^ "max * 

for pi 1. Here the exponent m is a constant and the 
function f(z) — > as z — > 0. Since the cost function 
attains its minimum at the boundary point we may set 
m = 1. Moreover, the latter term on the right-hand side 




FIG. 2: Characteristic form of the cost function, Eq. 12.91 . 

of Exp. (|2.7(l can be interpreted as a certain "interaction" 
potential between the following and lead cars which is 
long-distant one for m = 1 . A detailed analysis of effects 
caused by the value of the exponent m requires an indi- 
vidual study. Here the value m — I is actually chosen as 
just a simple reasonable assumption. Since the effect of 
the surrounding cars is depressed for small values of the 
car velocity v the function f(z) is also to attain its min- 
imum at z = 0. Therefore, f(z) = z 2 as z <C 1 should 
hold. In other words, inside a certain neighborhood of 
the origin {p = 0, v = 0} the following expansion 

r{p,v)=T{Q,v) + ^—pl (2.8) 

max 

can be adopted. 

Taking into account these speculations about the be- 
havior of the cost function !F(p, v) caused by variations 
of both its arguments, the following simple Ansatz will 
be used subsequently: 

\ Vmax / "max 

Figure El displays this function. It has only one global 
minimum at p = and v = i^max- For a fixed velocity 
v it attains a local minimum at the boundary p = 0. 
Ansatz ()2.9|l generalizes the adopted assumptions about 
the cost function. The former term takes into account 
relations l|2.4|l . I|2.5|l . the latter one is based on approx- 
imation H2.8fl . Of course, the parabolic Ansatz l|2.9|l is 
approximate only in the limit pi <C I . Since this is the 
main region of interest, this does not constrain its use- 
fulness. 

To deal with the car dynamics we should construct the 
cost functional £{h(t)} for the car motion paths {h(t)} 
treated now as continuous functions of time t. We note 
that the time dependence of headway distance hit) gives 
us the complete information about the car dynamics due 
to the relationship dh/dt = V — v. Leaping ahead, we 
say that facing this problem it is necessary to introduce 
additional notions. First, we should expand the phase 
space in describing the car motion state because transient 
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processes are now the subject of consideration. Second, 
drivers plan their behavior for a certain time in advance 
instead of simply reacting to the surrounding situation. 
So we should specify the region inside which a driver can 
monitor the traffic flow evolution and, thus, plan driv- 
ing her car. A priority functional similar to Eq. H2.2JI 
must span the time interval corresponding to this region. 
Beyond it the contribution of the path fragments to eval- 
uating the car motion quality at the given moment of 
time has to be fairly minor. 



III. RATIONAL DYNAMICS OF CAR MOTION 

A. Cost functional and the extremal equation 

Dealing with transient processes in the car motion we 
should consider once more the collection of phase vari- 
ables characterizing the cost of car motion at the current 
moment of time. Keeping in mind conventional driver 
experience, we will expand the current state of car mo- 
tion given by headway h and velocity v with the car 
acceleration a. This is essential because a driver can- 
not change the position and velocity of her car imme- 
diately, they vary continuously in time and contain no 
sharp jumps. Conversely, a driver controls the acceler- 
ation directly governing the car motion. Besides, she 
can change the acceleration practically without delay be- 
cause in the present analysis it is quite reasonable to ig- 
nore time scales related to physiological properties of the 
driver or to the mechanical properties of the car. So, the 
cost function T (h, v, a) for the motion state {h,v,a} 
can be written: 



F d {h,v,a) =J r (h,v) + 



§ 2 

max 



(3.1) 



where the time scale r > 1 s characterizes the acceler- 
ation capability of the car. In writing this expression 
we have assumed that driving without acceleration is 
preferable. Then, the cost function T& (h, v, a) had been 
expanded into a Taylor series with respect to a, keep- 
ing the leading term only. According to the result to 
be obtained the time scale r entering expression (|3.1() 
and the one used in equation (jl.4a.fl are practically the 
same. It should be noted that confining ourselves to ex- 
pansion (|3.1|) with respect to a the difference between 
acceleration and deceleration processes has been lost. In 
reality, they are different. Ignoring this difference leads 
to models where cars can crash. It is possible to take 
into account this effect using the approach under devel- 
opment, which, however, is worthy of individual investi- 
gation and will be done somewhere else. 

For a real driver, various thresholds in the driver recog- 
nition of hazards and obstacles exist [4l|. One of them 
is the distance A at which a driver can recognize the be- 
havior of other objects. This distance is usually related 
to the threshold of the visual angle 9 C subtended, for ex- 
ample, by vehicles ahead (Fig. [3}. The value of 9 C can 
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FIG. 3: Illustration of the recognition distance A and its re- 
lationship with the limit angle of perception 9 C . 

be estimated as 9 C ~ 15-30 min of arc ji^. The criti- 
cal angle 9 C , the characteristics height g of cars, and the 
corresponding mean distance A are related by 

A 

This allows for the estimation of the recognition distance 
as A ~ 200-400 m setting g ~ 2 m and 9 C ~ 30-15 min- 
utes of arc. As previously j^, we relate the driver an- 
ticipation with the region of size A in front that is clearly 
observable and where she can recognize the car behavior. 
A driver plans her motion based on the information re- 
ceived by monitoring traffic flow inside the observable re- 
gion. Under normal conditions this region should enable 
her to govern the motion effectively, for example, to de- 
celerate in advance avoiding a possible accident. There- 
fore its size A has to meet the inequality A > i5 max T. 
Leaping ahead, we introduce the value 



'max T 

a = — ; — < 1 



A 



(3.2) 



treated as a small parameter in the theory to be devel- 
oped. In particular, for r w 1 s, $ m ax ~ 100 km/h, and 
A « 300 m we have a rj 0.1 

The size A of the recognition distance can be estimated 
using another argumentation. If ~ 10 s is the typical 
deceleration time from the velocity i? max attained on the 
given empty road to zero value then the estimate A ~ 
i9 max Trf should hold. For ^ ma x ~ 100 km/h we have again 
A w 300 m. 

Now the cost functional of car motion can be written 
in an integral form as expression (|2.2f) containing the in- 
tegrand H3.1|) . The car is assumed to be located at point 
x along the road and to move with speed v at the current 
moment t of time. The possible paths of further motion 
{r(t, t), t > t} form the set &(t,x,v) on which the cost 
functional is defined. Here and below Gothic letters will 
be used to label the path variables. The function {r(t, t)} 
allows for a derivation of all the dynamical variables, the 
headway distance t)(t), the car velocity u(t), and the ac- 
celeration o(t) of the trial car path. The cost functional 
evaluating the quality of the path {r(t, t)} can be written 
as 



£{r} = 



poo 


r t(t)-xj 


L exp 


A 



p DO 









f d (\),u,a)dt, (3.3a) 



^ d ((],u,a)dt. (3.3b) 
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Forms (|3.3a|l and (|3.3b|) correspond to the driver pre- 
diction with spatial or temporal limitations, respectively. 
The cost functional of form Ij3.3a|l was used in [231 ■ Sub- 
sequently, functional (|3.3bfl is used because its form sim- 
plifies the mathematical manipulations when the velocity 
V of the lead car is fixed. In particular, dealing with func- 
tional (|3.3b|) we can specify directly the set of trial paths 
&{t, h, v) using solely the time dependence of {()(t, £)} of 
the headway distance. The path variables u(t, t), a(t,t) 
are completely determined by the dependence {[)(t, t)}: 



u(t,t) = V- 



dt 



a(t,t) = - 



d 2 t)(M) 
dt 2 



(3.4) 



Furthermore, the cost functional i|3.3bjl matches the 
driver prediction with temporal limitation, the mecha- 
nism which also allows for the interpretation of the recog- 
nition region size in terms of A ~ t? max T,;. However, both 
of the cost functionals lead practically to the same results 
for the analyzed situation. 

Each one of trial paths originates at time t and starts 
from (h, v) on the phase plane determined by the head- 
way distance and velocity of the car at the current time t. 
Since we investigate the car motion inside traffic flow but 
not processes of leaving it we assume that time variations 
h(t) of the headway distance are bounded. Therefore the 
trial paths fulfil 



KM) !*=* = &(*)> n(t,t)\ t=t = v(t) 



(3.5) 



and do exhibit bounded variations only as time goes to 
infinity. In other words, there is a constant C > with 



\\)(t,t)\<C for t>t. 



(3.6) 



In what follows, the driver behavior will be described 
by the extremals of the cost functional (|3.3b|) . Using the 
standard variational technique and taking into account 
conditions l|3.5|l . (|3.t)|) the governing equation for these 
extremals can be derived: 

dt 2 \dt a ' + X 2 

- j { d u F d + jd u F d - d^ d = . (3.7) 

To study the spectral properties of Eq. (|3.7|1 it is lin- 
earized around the stationary solution (hy, V). Its eigen- 
functions can be found by the Ansatz 



t 



f) C (t) cx exp -C- 



The obtained eigenvalue equation is given by 

(C + 0) 2 C 2 -A(C + ^)C+^ = O, (3.8) 
with the coefficients 

<$> = ^— a , A = itfL^F > , (3.9) 

"max ^ 

n = 2r 2 tf 2 nax (d 2 h T - jd h d v A (3.10) 



and the derivatives are taken at the stationary point 
(hy,V). The inequality d 2 T > is supposed to hold 
beforehand, in particular, it is the case for the cost 
function (12.91) . Equation l|3.8l) possesses four roots, one 
pair {£+, £_} of them have positive real parts, the other 
{Cj-jC— } have negative ones. Since the eigenfunctions 
with the eigenvalues {C+, C~} diverge as t — > oo we must 
omit them by virtue of condition 1|3.6J) . Roughly speak- 
ing, these divergent eigenfunctions describe the process 
of a driver leaving traffic flow, for example, to stop the 
car. Such processes are not under consideration. The 
former pair of eigenvalues are given by: 



c± 



2 2 



(3.11) 



So, in a small neighborhood of the stationary point 
(hv,V) of equation (|3.7|l any extremal f) pt(t,t) can be 
written as 



&opt(t)*) = h v + h+expl-C 



h- exp — ( 



t-t 

T 
t-t 
T 



, (3.12) 



where h + and are some constants. 

Equation l|3.7|l is of fourth order, so its general solution 
is specified by four conditions. Solutions that diverge as 
time goes to infinity have to be omitted (they are related 
to {C+)C-})- Their divergence is due to the exponen- 
tial cofactor in integral (|3.3b(l so it is retained beyond 
a small neighborhood of the stationary point (hy , V) . 
Therefore, only two conditions are needed to specify the 
desired solution of Eq. (|3.7|l . In particular, the initial 
headway distance h and the car velocity v determine it 
completely. The latter statement can also be proved for 
the more general form l|3.3a|l of the cost functional |33J . 



B. Rational driver behavior and Nash equilibrium 

In the previous section the dynamical cost functional 
was stated. At the next step the driver strategy based on 
this evaluation of the motion quality should be described 
(Fig. 2}. Following |33|. it is supposed that the driver, 
first, in planning the further motion chooses the path 
f)opt (t, t) minimizing the motion cost functional (|3.3b|) , 



driver: s-> f) op t(t, t) 



min £{t)(t,t)} . (3.13) 

ti(t,t)e&(t,h,v) 



To do this the driver "solves" equation (|3.7|l subject to 
conditions (|3.5|l . (|3.t)|) and, thus, get the optimal path 
of further driving, f) opt (t, t) . As noted in the previous 
section the optimal path choice is completely determined 
by the current car velocity v and the headway distance 
h. The terminal conditions, i.e., the goal of reaching the 
steady-state motion is implied. Since the driver controls 
the car motion through choosing the adequate value a 
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FIG. 4: Illustration of rational driver strategy. The left panel 
corresponds to the case where the driver's estimate of the 
car motion quality is not perfect. This is described by the 
cost function T(§,u,i \ t) containing the current time t as 
its argument. So the found optimal paths f) op t(t, t) depend 
on the time t when the driver plans the further motion. The 
right panel presents the ideal case where the driver is able 
to evaluate the motion cost precisely and the corresponding 
cost function J-"(fj,u, t) does not contain the time t. So, the 
optimal path family degenerates into one curve. It should be 
noted that the latter is correct even if the time t enters the 
cost function via variations of the lead car velocity. 



of acceleration she has to "find" the second derivative 
of f)opt(t, t) with respect to the former argument t (see 
Exp. (jsQJ) and, then, to "calculate" the result at the 
current time. Therefore, 



a(t) 



lim 

t-tt+0 



d 2 fro P t(t, t) 
dt 2 



(3.14) 



relates the current car acceleration a to the current head- 
way distance h and the car velocity v. 

Formula l|oM4fl describes the driver's choice at the cur- 
rent moment t of time. To convert it into the governing 
equation of car motion we adopt the second assumption 
that the driver performs this choice continuously: 



a = K{h,v,V) 



(3.15) 



All the car following models based on this equation, 
but with different particular forms of the cost function 
J rd ((),u, a) may be categorized as the rational car dy- 
namics approximation. 

Equation (|3.15l) holds even if the cost function 
T (h, u, t | t) depends explicitly on the time t at which the 
driver plans her further motion. Then, the driver evalu- 
ation of traffic flow state will change in time. Therefore, 
the resulting path {h(t)} of the real car motion envelops 
the family of the optimal paths {f)opt (t, t)} generated at 
different moments t of time (Fig. 0] left panel) . 

The present paper assumes driving to be perfect. The 
driver's choice is based on the precise knowledge of the 
cost function T d (h,v,a,t) depending solely on the cur- 
rent headway distance h, car velocity v, acceleration a, 
and, may be, the current time t. It does not depend 
on the moment of time when the driver evaluates the 



possible paths of her further motion. This assumption 
leads to what is called the Nash equilibrium in game the- 
ory [43. If, at time to the driver has found the opti- 
mal path f) op t(t, to) of the car motion (which meets equa- 
tion i|3.7|) and does not contain explicitly to), then any re- 
computation a certain time f > to later gives the same re- 
sult, t) opt (t, t) = f)opt(Mo) f° r t > f . In other words, if at 
time to the driver has chosen an optimal path f) op t(t, to) 
then the further motion will be described by it indepen- 
dently of either the driver follows it without correction 
or optimizes the car motion continuously (Fig. 21 right 
panel) . 

The Nash equilibrium in the driver strategy enables 
us to conclude that equation l|3.7|) describes the real car 
dynamics, not only the imaginary paths existing in the 
driver's mind during her planing of the further motion. 
In particular, it can be rewritten in a form containing the 
real acceleration o, velocity v, and headway distance h: 



d 2 a 



V da V 2 
dt 2 Jdt + l? a 

-^(jt^-jd^F + ih.lF) =0. CI hi) 



In the vicinity of the stationary point (hy, V) its solution 
is actually given by Ansatz (|3.12[) . which immediately 
leads us to the following expressions for the amplitudes 



h+ = 
h- = 



T(v-V)-(_(h-h V ) 

C+-C- 

( + {h~hy)-T{v-V) 

C+-C- 



(3.17a) 
(3.17b) 



and relates the car acceleration a to the car velocity v 
and the headway distance h: 



(C+ + C-) 



(v-V) 



CC+ (h-h v ) 



(3.18) 



(C+ + C-) 

These results are analyzed individually. 

1. Optimal driving condition 



The stationary point (h = hy, v = V, a = 0) of equa- 
tion (|3.16(l gives the headway hv which the driver chooses 
in order to follow the lead car at speed V. In particular, 
the expression 



d h T - -rd v J~ = 



(3.19) 



specifies the relationship between the values of the head- 
way h and the car velocity v for the stationary traffic 
flow imitated by the given car following problem. Solv- 
ing Eq. I|3.19fl for the car velocity v the optimal velocity 
approximation is obtained, v = $ pt(M- I n particular, 
for the specific form (|2.9|) of the cost function J-(h, v) we 
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FIG. 5: The kinetic coefficient g^ and the ratio t v /t vs f2. In 
drawing these curves Exp. 13. lit with A = 1 and = was 
used. The solid curve is 5^(0), the dashed curve is t v /t(Q). 



immediately get for I <C h <C A 

$opt(h) — $max ^2 i £)2 ' 

where the spatial scale I? is given by expression 
A/" 



D = \l — , ;<l><a. 



(3.20) 



(3.21) 



Relation (|3.2U|) or similar sigmoid functions are widely 
used in current literature. In particular, for l = lm and 
A = 300 m Exp. I|3.21|l gives the estimate D w 12 m 
typically ascribed to the spatial scale D. 



2. Linear governing equation for the car following problem 




1.0 1.0 2.0 3.0 4.0 5.0 

dimensionless time, t/(2r v ) 




0.0 1.0 2.0 3.0 4.0 5.0 

dimensionless time, t/(2r v ) 

FIG. 6: The relaxation curves of car dynamics for different 
values of the parameter f2 (labelled by numbers in squares). 
The solid curves display the relaxation caused by a pertur- 
bation in the headway, the dashed one shows the same for a 
velocity perturbation. The variations in the headway distance 
and relative velocity are measured either in units of the initial 
values Sho and Svo or 2t v 5vq and 5/io/(2t„), respectively. 



By introducing the time scale t v = t/(£ + + C— ) an( i 
the coefficient 



9h 



C-C4 



(C+ + C-Y 

the car dynamics equation (|3.18(l can be rewritten as 



(3.22) 



1 



(v-V)- g h 



(h-h v ) 



(3.23) 



Equation (|3.23(l plays a significant role in models that 
describe non-rational driver behavior |44j , where expres- 
sion (|3.23l) specifies an optimal acceleration that could 
be chosen by rational drivers and the parameter gh was 
introduced phenomenologically. Here, gh can be com- 
puted using the cost function l|2.9|) . It is plotted as well 
as the ratio t v /t as a function of the variable f2 (Fig.EJ). 
This representation is due to the fact that the coefficient 
A — 1 (for the cost function (|2.9() A ~ 1 for I < h < A), 
the coefficient <C 1 by virtue of the adopted assump- 
tion (|3.2|1 . whereas the coefficient SI varies around unity. 
Below, the latter will be justified and the dependence of 
the coefficient Jl on the car motion state will be analyzed. 
As seen in Fig.GDthe time scales t v and r practically co- 
incide with each other and the parameter gh is a small 
value. 



3. Relaxation curves vs system parameters 

Using again the cost function Eq. I|2.9|) , and taking into 
account I <C hy <S A and A ~ 1, the optimal velocity de- 
pendence ^opt(^v) is given by Eq. (|3.20(l . The parameter 
f2 depends as follows on the car motion state 



n = Aa 



Td-d opt (h) 



dh 



= 8a (, 

h=h v {n 



XD 2 h v 



D 



2\2 



It attains the maximum 
3%/3(rtf n 



2DX 



3V3a 2 A 
2 D 



(3.24) 



(3.25) 



at h n = D/x/3. Using t ~ 1 s, tf max ~ 100 km/h, 
D ~ 15 m, and A ~ 300 m, respectively, tt max ~ 0.5 is 
obtained. Of course, this is not a precise numerical value, 
but in general f£ ma x ^ 1 may be assumed. 

The precise value of f2 max matters. According to 
Sec. lIII^l the car dynamics is characterized by two eigen- 
functions with eigenvalues given by Eq. H3.11fl . When 
ft < 1, the relaxation is a pure fading process character- 
ized by time scales t/£ + and r/C-- Otherwise, Q > 1, the 
relaxation process is characterized by the complex eigen- 
values, However in this region their dependence 
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on O is sufficiently weak, |C±| ~ f° r ^ ^ 1 as follows 
from Exp. (|3 . 1 1|1 . Therefore the car relaxation processes 
should be rather insensitive to the value of fl as it varies 
from unity to about ten. Figure El demonstrates this fact 
showing the relaxation curves induced by deviations from 
the stationary values of the headway and the car velocity 
individually. All the corresponding curves lie near each 
other, only for O = 0.1 a larger deviation is visible. This 
is because for such values of f2 the eigenvalues t/(± dif- 
fer essentially. By contrast, the cost of acceleration, i.e., 
the cost of correcting the car motion depends strongly 
on f2 because this component of the cost function (|3.1|) 
varies linearly with Q to the first approximation. Since 
the maximum r2 max of the parameter f2 is attained at 
h ~ D this means that the greater the value of O, the 
larger is the cost of the car governing with respect to the 
cost of motion in its own accord. 
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FIG. 7: The weight coefficient n of the combined model ll.4al 
vs the stationary headway hv- For comparison, the dashed 
curve shows the optimal velocity, Eq. 13.2011 . 



4- Following-the-leader model vs optimal velocity model 



Hypothesis about the value of "d ma x 

Keeping in mind the aforementioned speculations let 
us assume that the composed parameter f2 max takes a 
certain fixed numerical value, f2 max = f2 c < 1. It leads 
immediately to a certain relationship between the initial 
quantities of the model, i9 max , A, r, I. In this collec- 
tion only the optimal velocity i? max on empty road is 
determined solely by the driver behavior. The other pa- 
rameters, namely, the size A of recognition region, the 
characteristic time scale r of the car ability to acceler- 
ate or decelerate, and the characteristic headway I that 
drivers can control reliably and should be about the typ- 
ical headway distance in dense jams are determined by 
other mechanisms, e.g., car properties, driver physiology, 
visual conditions and the like. At first glance such re- 
lations seem to be impossible, because a driver chooses 
the optimal velocity on empty road depending on the 
road conditions, speed regulations, weather, etc. Nev- 
ertheless, there is an additional factor affecting the car 
motion on empty roads, the driver experience. If the 
driver preference for the empty road velocity stems from 
here experience of driving in low density traffic flow then 
such relations can exist. 

In the given model this assumptions lead to the ex- 
pression 



As discussed already there are two stimuli for a driver 
to change the current motion state. One is the velocity 
difference v — V between the leading car and her car. 
The second is the difference between the current speed 
v and the optimal speed $ pt(^) as function of headway. 
Models like (|1.4|) take into account both of them, leaving 
the determination of the weight coefficients to appropri- 
ate empirical or experimental data. The results obtained 
so far actually allow to predict these coefficients. 

This can be demonstrated with the combined 
model (|1.4a|) . Linearizing Eq. (|1.4a|l around the sta- 
tionary point (hy,V) and comparing the result with 
Eq. ijXTgj) we get 



C-C4 



(C+ + C- 



opt 



dh 



h — h\ 



(3.27) 



This expression gives k as function of the car motion 
state. Figure [7| plots this dependence for the cost func- 
tion (|2.9|l . As seen in Fig. [7| the weight coefficient k is 
small for all interesting values of the headway distance 
in the car following regime. Only for small headways 
hy ~ I (dense jams) or for large headways hy exceeding 
D substantially (free flow) k approaches unity. 



5. Different types of the car dynamics 



L'max 



,1 \ 3 



(3.26) 



Naturally, the specific form of relation (j3.26|l will change 
if, for example, another Ansatz of the cost function is 
used, however, its qualitative features should be retained. 
These speculations, definitely, are not enough to regard 
relations similar to formula l|3.26|l as an established the- 
oretical result. We have presented it for discussion only. 



As the headway distance hy varies in the interval I <C 
hy <C A the coefficient Q(hy) changes essentially, leading 
to qualitatively different car dynamics. As mentioned 
before, when Q(hy) > 1 the relaxation exhibits damped 
oscillations around the stationary point (hy, V). This is 
due to the eigenvalues {C+,C-} ( sce Exp. (|3.11|l ) having 
non-zero imaginary parts. Since for Q(hy) > 1 their real 
and imaginary parts are of the same order the car motion 
relaxation is characterized by one time scale about one 
(in units of r). This is true also for Q(hy) < 1 because 
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FIG. 8: The ratio of time scales r = IC-I/IC+I °f the corre- 
sponding eigenfunctions describing the car motion relaxation 
as a function of SI or the headway hv when il c is fixed. Com- 
putations in this figure used cost function 12.91 and = 0. 



C + « £_ there, although the dynamics is a pure fading 
process now. For still smaller values of tt(hy) ~ 0.6 the 
ratio r — |C-|/|C+| of these eigenvalues becomes smaller 
than one half, which may be defined as the two-scale 
regime of the dynamics, see Fig. |H] 

The change between the two-scale and the one-scale 
regime corresponds to a value of hy ~ 1-1.5 D. Below 
this value, the dynamical behavior is one-scale, above 
it is two-scale and will be called "fast-and-slow" in the 
following. 

Concerning with the "fast-and-slow" dynamics the ve- 
locity relaxation and the headway relaxation can be an- 
alyzed individually. According to Exps (|3.17|) the initial 
difference h — hy contributes mainly to the eigenfunction 
with the eigenvalue Thereby, the velocity difference 
v — V contributes mostly to the amplitude h+. The am- 
plitude /i_ also contains the term of the same magnitude, 
however, the time scale t/£_ on which the corresponding 
eigenfunction varies is much bigger than the time scale 
t/( + of the other eigenfunction. So, the velocity relax- 
ation falls on the first eigenfunction. Thus, the velocity 
difference v—V disappears practically completely during 
the time r/£ + , which forms the "fast" stage of the car re- 
laxation. At the next "slow" stage of duration t/£ + the 
headway deviation from the equilibrium value hy disap- 
pears. To summarize, the "fast-and-slow" car dynamics 
is a two-stage process where the velocity difference be- 
tween the cars is eliminated first. Later on, the headway 
is optimized. While this is being done, the resulting ve- 
locity difference is not essential and the driver can govern 
the car motion without the necessity of responding fast. 

Another important characteristics, separating dense 
traffic and quasi-free flow can be derived as follows. The 
solution of the eigenvalue equation i|3.8|) depend actu- 
ally on two values, 4> and f2, since A ~ 1 for the cost 
functional (|2.9(l . The parameter <fi < a <C 1 by virtue 
of the adopted assumption l|3.2[l . The maximum of f2 
attained at hy — £)/v3 is much larger than a 1 as it 
results from Exp. (|3.25|) . However, as the headway dis- 
tance increases the value Vl(hy) decreases as h v 3 (see 
Exp. H3.24[l ) whereas <f> — > a. So there is a value h c of the 
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FIG. 9: Possible types of the car motion state classified us- 
ing different properties of the car dynamics and their mutual 
arrangement depending on hy- 



headway distance at which both these terms contribute 
to the eigenvalues to the same extent. To find h c consider 
h ^> D where ft(hy) <C 1 and <\> w a, thereby, £ + ~ 1 
and 



C- 



2(<j+ y/a 2 + n(h v )) 
by virtue of formula H3.11|) . So, h c is specified by 



(3.28) 



n(h c 



i 

//,. =21 -^y ; />. 



(3.29) 



The value h c divides the headway distance region into 
two parts, see Fig. When hy > h c the velocity of 
the following car is close to $ max , so this type of car 
motion will be referred to as the quasi-free flow. In this 
case the eigenvalue equation (|3.8|) cannot be simplified, so 
to describe the quasi-free flow the cost functional l|3.3b|) 
only in its full form may be used. In the opposite case, 
h <C h c , which will be called dense traffic mode, the term 
4> is ignorable, reducing the eigenvalue equation (|3.8|) to 



C 4 -AC 2 



1 



= 0. 



(3.30) 



Moreover, in the latter case the description of car dy- 
namics can be reduced to the standard form of classical 
mechanics, enabling us to analyze the nonlinear stage of 
the car dynamics. 



IV. NONLINEAR CAR DYNAMICS 

A. Effective cost functional for dense traffic flow 

In the dense traffic limit h <C h c , the variational prin- 
ciple based on optimizing the cost functional l|3.3bj) can 
be simplified essentially. In the given limit the char- 
acteristic time scales of the car dynamics are r and 
t / ^jQ,(hy), with both of them being small in compari- 
son with A/$ max - The latter follows from the adopted 



assumption (|3.2|l and the condition Q(hy ) 
h < h c - 



> a' 



for 



11 



In this shown in Appendix [5] the cost func- 

tional (|3.3b|) can be replaced by the following effective 
functional whose integrand does not contain a time de- 
pendent factor 



C{h(t)} 



T e s (h, v, a | V) dt' . 



(4.1) 



where the integrand, which will be called the Lagrangian 
of the car dynamics is given by 



T cS (h,v,a\ V) = T d (h, v, a) 

'V 



h + v) d v !F\ 



V,hv 



(4.2) 



Interestingly, T e s (h, v,a\V) contains the lead car veloc- 
ity V as a parameter and attains its extremal value with 
respect to h at the point hy corresponding to the opti- 
mal driving with the velocity V. The term vd v !F\v,h v 
has been introduced for the sake of convenience only, it 
does not affect the extremal equation but enables the 
Lagrangian to attain a minimum with respect to the car 
velocity v at the stationary point (/iy, V). 

Both functional l|3.3b|) and l|4.1|l possess the same 
extremals to the first order in the small parameter 
In 



a I y/Q(hv)- In particular, the 
tional 1)4. 1|) meet the equation 



extremals of func- 



2t 2 d 2 a 

^max dt 



0,, 



V 



h\=0. 



(4.3) 

It corresponds directly to the initial full equation (|3.1()|) 
where the second and third terms in the former paren- 
theses are ignored whereas the second term in the latter 
parentheses is replaced by its value taken at the station- 
ary point. It is justified because these terms are due to 
variations of the time dependent cofactor. 

Keeping in mind its following applications we rewrite 
Lagrangian T c s (h, v,a \ V) also in the form 



^■ eff (h, v, a | V) = 



$2 
max 



+ F {v) + T- mt (h,v | V) . (4.4) 



In particular, for the cost function (|2.9|l we get 



H{v | V) 
JWM I V) 



(v - Vf 
d 2 

max 

nh 2 v 

max 



l - 



v 2 

d 2 

max 



v 2 h v 



h 

hv 



(4.5) 
(4.6) 



In obtaining formula l|4.6|) we have taken into account 
the expression (|3.2U|I relating the velocity V to the head- 
way distance hv, expression (|3.24() . and omitted some 
insignificant terms. The latter term in expression (|4.5|) 
can be also omitted because it has no effect on the ex- 
tremal equation. 

The most essential feature of the effective cost func- 
tional H4.1fl is the absence of a time dependent cofactor. 
This enables us to reformulate the car dynamics in terms 
of autonomous Hamiltonian equations. 



B. Hamiltonian description of car dynamics 

Finding the extremals of functional 14.1fl can be done 
as follows. The phase space {h 7 v,a} is expanded to 
{h, h, v, v} by adding the relationship between the head- 
way distance h and the car velocity v as an additional 
constraint: 



minimize 



Jo 



T e s(h,v,v)dt 



subject to the dynamical equation 



h = V 



and the initial and final conditions, 



h(0) = h , 
h(oo) = hy , 



v(0) = v , 
v(oo) = V . 



(4.7) 



(4.8) 



(4.9) 
(4.10) 



For the sake of simplicity the parameter V has been omit- 
ted in the list of variables of J- e s, the current time is set to 
zero. By using a Lagrange multiplier, the problem 14.7(1 
is rewritten in the standard form. Namely, the extremals 
of 



T e tt (h, v, v) + p (h + v - V) 



dt 



(4.11) 



are sought. They are defined on the extended set of 
variables {h(t), v(t),p(t)} and subject to conditions (|4.9|1 . 
(|4.1U|) . The extremals of functional l|4.11|l obey the clas- 
sical Lagrange equations. 

In constructing the Hamiltonian TL{h, v,p, q) that pro- 
duces the same extremal equations the Pontryagin tech- 
nique is used 0| . Introducing the new variable q 



dT cS (h,v,v) 



dv 

and solving for v, i.e., finding v 
Hamiltonian is written as 



(4.12) 



v(h,v,q), the desired 



H(h, v,p, q) = qv - T cS {h, v, v) - p(v - V) . (4.13) 

It can be demonstrated directly that the desired ex- 
tremals meet the standard Hamiltonian system of equa- 
tions 



h 
P 



~d h n. 



d q H, 

-d v n. 



(4.14) 
(4.15) 



Again, the optimization of functional l(4.7|l actually 
leads to a boundary value problem because the extremal 
is specified by both the initial and final conditions 1)4. 9|) . 
H4.10|) . The Hamiltonian approach H4.14f) - I|4.15|l also 
shares this property. Namely, the initial values of the 
variables h and v are given by conditions (|4.9|l . whereas 
the initial values po and qo of the quasi-momenta p and 
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q should be chosen such that the system tend to the sta- 
tionary point (hy, V) as t — > oo. It is the same situation 
that we met in dealing with equation H3.16JI 

However, the Hamiltonian description has a certain ad- 
vantage. First, the Hamiltonian itself is conserved, 

dM(M,P,g) =0- (416) 

at 

Second, an additional autonomous first integral of the 
system (|4.14|) . (I4.15|l . can be found, at least, for a certain 
stage of the car dynamics. Fortunately, this is the case 
for the "fast-and-slow" car dynamics analyzed below. 

Note, that this Hamiltonian description of the car mo- 
tion relaxation towards the stationary state does not con- 
tradict the conservation of phase volume being the gen- 
eral property of Hamiltonian systems. The matter is that 
the relaxation process is described by only one path lead- 
ing to the stationary point of the saddle type. Other 
possible paths are not considered. 

Physical meaning of the Hamiltonian variables 

The two Eqns (|4.14|) can be understood readily. The 
equation for h is just h = V — v, while the equation for 
q, or, what is the same, expression <|4.12[1 shows that the 
quasi-momcntum q is proportional to the car acceleration 
a = v: 

2r 2 

q= w a. (4.17) 

The physical meaning of the quasi-momentum p is 
more complex. To clarify it, rewrite the equation 
line H4.15fl m terms of partial derivatives of the La- 
grangian T c s{h, v, v). Namely, the definitions of the 
quasi-momentum q <|4.12[1 and the Hamiltonian 14.13[) en- 
ables us to represent these equations as 

q = d v F eS +p, (4.18) 
P = d h T cS . (4.19) 

In particular, expressions l|4.17[l - i|4.19|l immediately lead, 
as it must, to equation (|4.3H . 

Expression l|4.18|l demonstrates that the quasi- 
momcntum p includes the rate d of acceleration varia- 
tions. In this way the variable a implicitly enters the 
Hamiltonian 7i(h, v,p, q) specifying a hyper-surface in 
the four-dimensional phase space {h, v, a, a}. Even more, 
by ignoring the dependence of !F c s(h, v, v) on the head- 
way distance h, the velocity relaxation towards the sta- 
tionary value V is governed by the equation 

q = d v T cS ■ 

It results from the effective cost functional (|4.1|l where 
the car velocity v is treated as a primary argument. 
As will be seen below, there is a stage of the car dy- 
namics where this assumption is justified and the quasi- 
momentum p = is conserved. Only the dependence of 



the Lagrangian T c s(h, v, v) on the headway distance h, 
see Eq. (|4.19|) . causes variations in the quasi-momentum 
p. Since the velocity control is of primary importance in 
the driving strategy, the variable p can be regarded as a 
certain measure of the necessity to control also the head- 
way distance when eliminating the velocity difference. 

C. "Fast-and-slow" dynamics 

For £1 <C 1 the last term on the right-hand side of 
Exp. I|4.4|) is small. Since the headway h enters the 
Hamiltonian l|4.13[l exactly via this term time variations 
in the quasi-momentum p are retarded by virtue of the 
first equation at line l|4.15[) . The same follows from 
Eq. H4.19fl . The other variables h, v, q can vary substan- 
tially on scales independent of the value Q -C 1. Then, 
the car dynamics can exhibit multi-scale relaxation. In 
particular, Sec. IIII B 51 demonstrated this fact, namely, 
it has been shown that in the given limit the system 
dynamics comprises the fast and slow stages. Exactly 
the time scale of the latter stages depends on the value 
Q -C 1. During the former stage the quasi-momentum p 
should be practically constant, i.e., it is approximately a 
first integral of the system (|4.14|) . I|4.15[) . This property, 
together with the conservation of Hamiltonian (|4.13|) , al- 
lows the integration of the system of equations 14.14|l . 
(|4.15|l . During the following slow stage a quasi-stationary 
approximation can be used to analyze the system evolu- 
tion. 

The results below are exemplified with the cost func- 
tion (|2.9|l . They can be easily generalized to other cost 
functions. 



1. Fast stage 

At the zeroth approximation in the small parameter 
O the quasi-momentum p is a constant. During the fast 
stage mainly the velocity deference v — V is eliminated, 
so Eq. (|4.18|) formally describes the relaxation process 
of the car velocity v to the stationary value V. There- 
fore, by virtue of Exp. I|4.5|l . the quasi-momentum p takes 
zero value, p — 0. In this limit the Lagrangian compo- 
nent Fint(h, v | V) as well as the additive constant of the 
component Tq(v \ V) can be omitted. So, the Hamilto- 
nian (|4.13l) becomes 

H(v, a) = TL a {a) - H v [v) , 

where 

r 2 a 2 (v - V) 2 

H a (a) = — , H v (v) = [ — -L , (4.20) 

max max 

and the non-canonical variables {h, v, a} are used. 

Since the Hamiltonian H.(v,a) is conserved and the 
fast stage describes the velocity relaxation to V the car 
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collisionless j 
region I 
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initial car velocity 
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is worth of an individual consideration. Here we only 
touch on this problem by assuming the collision to hap- 
pen when the value /i* given by Exp. ((4.23|) becomes 
equal to zero. This assumption is justified as a rough 
approximation due to estimate (|4.24(l . see Fig. ^3 The 
boundary of the collisionless region can be shifted to the 
right substantially by dropping the assumed quadratic 
dependence of the cost function on acceleration. This 
will be done elsewhere. 



FIG. 10: Phase region of the initial values of the headway 
distance /io and the car velocity vo where the car dynamics is 
collisionfree (based on equation 14.231 '). 



dynamics obeys the equation 



H a (a) = TL v {v) . 



(4.21) 



This immediately gives the relationship between the ac- 
celeration a and the velocity v 



a= —(v-V) 

T 



(4.22) 



The sign in the latter equation has been chosen so to 
allow for the system relaxation. 

Equation 14.21(1 is, in fact, of the general form and the 
linear form of Eq. ((4.22(1 is due to the adopted quadratic 
Ansatz of the cost function but not a consequence of 
linearization. Besides, equality 1(4.21(1 can be read as fol- 
lows. The comparison of the Hamiltonian parts Ti a (a) 
and H. v (v) with the corresponding components of La- 
grangian l|4.4|l demonstrates that the function Ti a (a) ac- 
tually measures the cost of the car acceleration and the 
function Ti v (v) does the same with respect to the car mo- 
tion relative to the optimal driving conditions. Thereby, 
during the fast stage the driver corrects the car dynam- 
ics so that the cost of acceleration be equal to the cost of 
current motion measured relative to the stationary con- 
ditions. 

The headway distance h does not enter explicitly the 
governing equation of the fast stage. However, it varies 
during the fast stage and finally attains a value /i* dif- 
fering from the initial value Hq. In the adopted simple 
approximation of the cost function Eq. (|4.22|) enables to 
estimate easily the value h*. Namely, the direct integra- 
tion of Eq. (|4.22|) yields the relationship 



K = h Q - t(v q - V) . 



(4.23) 



Applying to formula (|4.33|) it can be seen that expres- 
sion H4.23fl holds until the headway distance h becomes 
too small and it is impossible to ignore the effect of the 
term J-- m t(h, v), i.e., when 



(4.24) 



The car dynamics in the region of small values of the 
headway distance when the probability of collision is high 



2. Slow stage 

When the system attains a quasi-equilibrium with re- 
spect to the car velocity v its further evolution is due to 
the direct dependence of the Hamiltonian 7i{h,v, q, p) on 
the headway distance h. It enters the Hamiltonian via 
!Fi n t(h,v | V). As a result, the velocity difference v — V 
and the acceleration a should be small. Keeping this in 
mind, Eq. ((4.18(1 can be solved for p. Substitution of 
the obtained expression into Eq. 1(4.19(1 . and the further 
linearization with respect to the variable v — V and the 
component .Fint results in 



ad 2 v T \ 



v=V 



(4.25) 



If t s is the characteristic time scale of the slow stage then 

. Tr , h — hy h — hy . h — hy 
(v - V) ~ , a— ^ — , gaa~ ; — . 

These estimates together with Eq. 1(4.2511 lead to 

t s - fT 1/2 T. (4.26) 

In particular, the rate q of time variations in the quasi- 
momentum q scales as q cx £1 2 for f2 — > 0. Since the 
Lagrangian J- C g has a minimum at h = hy the quasi- 
momentum p can be estimated as 

p ~ TgdhTint oc f2 1/2 (/i - h v ) 

by virtue of Eq. 14.19(1 . In the limit VI <C 1 the term q 
can be ignored in comparison with the quasi-momentum 
p. Then, Eq. ((4.18|) immediately leads to the relation 



P 



-d v T t 



eff • 



(4.27) 



which is no more than the standard expression for the 
momentum of a system described by the Lagrangian 
T c s{h, v, 0) with v = V — h. 

By the same reasons all the terms in Hamilto- 
nian 1(4.13(1 containing the quasi-momentum q may be 
omitted. Then the substitution of 1(4.27(1 into ((4.13(1 
yields the Hamiltonian of the slow stage 



H(h,v)=H v (v)-H h {h), 



(4.28) 



where 



H h (h) - F int (/i, V\V)- F m {hv,V | V) 
{h-h v f 



nh 2 v 



4i? 2 
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(4.29) 
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In the Lagrangian component T- m t (h, v \ V) the car ve- 
locity v has been replaced by its stationary value V, and 
some negligible terms have been omitted. Then the con- 
servation of the Hamiltonian TL{h, v) during the system 
dynamics gets the form 

H v (v)=H h {h) (4.30) 
which leads to the following governing equation 



v -V = ^§{h-hy)^. (4.31) 

The conservation law (|4.3UI) can be understood in analogy 
to that for the fast stage. However, during the slow stage 
of the car dynamics the car velocity plays the role of 
the control parameter. The driver changes the speed to 
correct the headway distance. Again the cost of deviation 
of the car velocity from the stationary value is chosen to 
be equal to the cost of the car motion measured relative 
to the optimal conditions. 

It should be pointed out that the slow stage of 
the car dynamics is governed by the conservation law, 
Exp. i|4.30|l . which contains only the headway distance h 
and the velocity v. At the first approximation, the ac- 
celeration a does not enter at all. In this meaning, the 
slow stage is similar to other physical systems. However, 
the stationary point of the car dynamics is a saddle point 
rather than a minimum. 

Up to now, the two different stages of car-following dis- 
cussed already in Sec. IIII B 51 have been identified with 
two different conservation laws. Both of them can be 
unified into just another effective conservation law inter- 
polating Exps. (l4~?Tj) and (|Q1 : 



equations for the following car motion the driver prefer- 
ence has been used to construct a cost functional. Its ex- 
tremals specify the optimal paths of the further motion 
for the following car. Applying to the general proper- 
ties of the driver behavior we analyzed the basic proper- 
ties of the cost function and proposed a simple parabolic 
Ansatz, which, nevertheless, catches typical features of 
traffic properties. 



The concept of a rational driver is formulated. It com- 
prises the assumptions that the driver follows the optimal 
paths and corrects the car motion continuously. In this 
case the optimal path is a Nash equilibrium of the sys- 
tem, which is defined as follows. If the driver has chosen 
an optimal path at a certain moment of time then no fur- 
ther correction of the car motion is necessary because it 
leads to the same result. As the consequence of the Nash 
equilibrium the extremals of the cost functional specify 
the real dynamics of cars with rational drivers although 
originally they determine only the imaginary paths in the 
driver's mind when planning the further motion. In this 
way we obtained several results. The optimal velocity 
approximation has been derived. The weight coefficient 
entering a car following model combining the "following- 
the-leader" model and the optimal velocity model has 
been found depending on the headway distance. As an 
important result it has been shown that the car dynamics 
can be categorized under different types according to its 
properties. First, we have shown that there can be two 
types of the car relaxation towards the stationary mo- 
tion, the mono-scale dynamics and the "fast-and-slow" 
dynamics. Second, we singled out the quasi-free motion 
and the dense traffic mode. 



H v (v)=H a (a)+H h (h) 



(4.32) 



Within the same accuracy it is possible to interpolate 
directly Eqs. (|4~22|) and JOTJ, leading to 



1 



a = 

T 



(v-V) 



~2T\~h 



(h - h v ) 



(4.33) 



Besides, the proposed interpretations of the conserva- 
tion laws 14.21|) and l|4.30|l enable us to formulate a gen- 
eralized principle of adequate control. It declares that 
the effective cost of correcting the car motion via chang- 
ing the car acceleration (fast stage) or the car velocity 
(slow stage) is equal to the cost of the current state of 
car motion measured with respect to the optimal driving 
conditions. 



V. CONCLUSION AND DISCUSSION 

A variational approach to the description of car dy- 
namics has been developed in case of a car following a 
lead car moving at a constant speed. To derive governing 



The variational technique for the latter mode can be 
simplified essentially. Namely, it is possible to reformu- 
late the cost functional so that it does not contain a time 
dependent cofactor. As a result, the autonomous Hamil- 
tonian description for the car dynamics has been con- 
structed. For the dense traffic mode the "fast-and-slow" 
dynamics has been analyzed for the nonlinear stage. In 
particular, different conservation laws for this stages have 
been found. A generalized principle of adequate car mo- 
tion control has been proposed. 



We assumed that the lead car moves at a fixed speed. 
If it is not the case two different situations should be 
singled out. For perfect drivers who can predict rigor- 
ously the motion of the car ahead the main results still 
hold. Otherwise the cost function will depend not only 
on the current time but also, what is crucial, on the time 
when the driver has started to evaluate the further car 
dynamics. This disturbs the Nash equilibrium and the 
car correction should be carried out continuously. Briefly 
this problem was studied in [33[ but it actually requires 
a more detailed investigation. 
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Beyond the rationality 

Real drivers have certain limitations. They are not ca- 
pable of finding the optimal path precisely and they can- 
not correct the car motion continuously. So, the concept 
of rational driver behavior is just the first approxima- 
tion of the real situation and deviations from this perfect 
behavior should be analyzed consistently. A first step 
towards this problem can be found in Ref. 0]. Here we 
justify some assumptions adopted there and substanti- 
ated them in a different way. 

Let us again concern the driver evaluation of the car 
motion quality. For a given path of the car's further 
motion {f)(t, t)} the cost functional (|3.3b(l is written as 



die 



r 2 a 2 
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d 2 . 



(5.1) 



Expression (|5.1|) can be represented also in the form 



£{*>} = 



dt 



■d 2 

max 
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max 



(h 
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- hy) + Cq 
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IS 2 



■d 2 



hy (hy 



(5.2) 



where S(i, t) describes the deviation of the given path 
from the stationary car motion trajectory, 

l)(t,t) = h v + S(i,t). 

Here, Cq is the cost of stationary car motion, 

V 



0„ 



1 



(5.3) 



In addition, when deriving expression (|5.1|l we have used 

I _ I IS IS 2 

Ti ~ h v h\ + hy {h v + 5) ' 

integrated various fragments of Eq. (|5.1|l by parts, as- 
suming I <C f) <CA, and dropped terms like where 
ever possible. 

When the driver plan her further motion the current 
values of the headway h and the velocity v are regarded 
as the initial conditions. So, when choosing the opti- 
mal path she can minimize only the first term £ c {f)} 
of expression Q5.2[) . This optimization is implemented 
through the adequate control over the car acceleration a, 
so, exactly the acceleration plays the role of the control 
parameter available for the driver actions. 

A real driver can only approximately evaluate the qual- 
ity of motion. Let us describe the threshold in the driver 
perception of the motion quality by 

Ahr = (5.4) 



where e c is a small constant, e c <C 1, because definitely a 
driver can recognize the difference in the state of staying 
and freely motion on an empty road. 

When the controllable part C c {()} of the cost motion 
is much smaller then the threshold (|5.4|l the driver actu- 
ally has no information of how to govern the car motion. 
In this case it is natural to assume that he will not do 
anything with respect to the car driving and, so, will fix 
the car motion at the current state, including the current 
acceleration a. Thereby the inequality 



£ c {f)} <C 



thr 



(5.5) 



determine the region in the phase space {h,v,a} inside 
which the driver cannot control the car motion. In eval- 
uating such a driver behavior with expression l|5.2[) we 
can formally treat 5, u, and a as constants independent 
of one another. So setting S = h — h v , U = v, and a = a 
and taking into account the relation between hy and V 
we get from condition l|5.5|l the approximate boundary of 
this region in the following form 



2 2 

r a 



+ (v-V) 2 + 



tt(hy) hy 
~^ 2 h 



(h-hy) 2 <e 2 c dl 



(5.6) 



We recall that the parameter entering Eq. I|3.23|l co- 
incides with v / f^/(2r). So we reproduced here the ex- 
pression for the rational driving boundary as has been 
introduced in paper 44] by a different line of reasonings. 
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APPENDIX A: ADMISSIBLE FAMILY OF THE 
COST FUNCTIONS 

A set S of one-dimensional paths {a:(0}lt*lo °^ car mo " 
tion is considered. All the paths are supposed to originate 
from one point xq at the initial time t = 0, i.e. x(0) = xq, 
and to go to infinity, x(t) — > oo, as t — > oo such that the 
car velocity v(t) — dx/dt is bounded during the motion, 
\v(t)\ < D (here D is some constant). 

It is assumed that the set § is ordered by a certain 
preference relation ^ and there is a functional £{x(t)} 
of the form 

/•OO 

£{x(t)}= dtw(x,t)J 7 (x,v,t) , (Al) 
Jo 

measuring numerically the preference relation, i.e. meet- 
ing the condition 

xi(t) ^ x 2 (t) <^ £{a:i(t)} > £{x 2 (t)} . (A2) 
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Here, first, the integrand F(x, v, t) depends on the car 
position on the road x(t), the velocity v(t), and, may be, 
the current moment t of time. It plays the role of cost 
function, i.e. evaluates the quality of the car motion state 
at current time. Second, the weight factor takes one of 
the forms 

w(x, t) = w a (x) > or w(x, t) = Wb(t) > , 

where w a (x) — > and Wb{t) — » as x — > oo or t — > oo. 
The functions w a {x) and Wb(t) describe the spatial or 
temporal limitations in the driver prediction, respec- 
tively. The examples 

w a (x) = exp and w b (t) = exp > 

where A and V are some constants, enable us to relate 
the cost functional IjAljl to ones used in the main text. 
The specific form of w a (x) and Wb(t), however, does not 
affect the results to be obtained in this Appendix. 

Let £{x(t)} be another cost functional of form (|A1(I 
with its own integrand !F(x,v,t) that also describes the 
same preference relation. The purpose of the present 
Appendix is to derive the relationship between the cost 
functions T(x, v, t) and !F(x, v, t). 

For this purpose the path set S is divided into classes 
{&} of equivalent paths, i.e. x\ (t) , x 2 (t) & 6 if x\{t) ~ 
x-z{t) or, what is the same, C{x\{t)} — £{x2(t)}. The 
cost functional £{x(t)} based on the integrand T{x, v, t) 
has to preserve this class partition. So the condition 
SC/Sx = must hold for any infinitesimal perturbation 
Sx(t) of a path x(t) within the class S[x(t)]. Using the 
Lagrange multiplier technique and taking into account 
the initial condition x(0) = xo (thus, Sx(Q) = 0) the 
latter requirement is reduced to the equation 

(A — vd v A) d x w — (d v A)d t w 

+ (d x A-j t d v Ajw=0, (A3) 

where the function A(x,v,t) := !F(x,v,t) — xT(x,v,t) 
and x is a certain constant for all the points (x, v = 
dx/dt) on the phase plane {x,v} related by the given 
curve x(t). Let us show that the value of h is constant 
for all the points of the phase plane. The cost functions 
T{x, v, t), J-(x, v, t) are chosen, thereby, if the quasi-local 
equation 1A3(I holds, the value of x must be determined 
by any fragment of the path x(t) under consideration, 
for example, by a fragment x*{t) := x{t)\^ G x(t)\^. 
For t < t\ and t > t% any point (x, v) can be joined to 
the fragment x*(t) by a certain path x(t) which must be 
characterized by the same value of k. The latter proves 
this statement because the time moments t\ < ti have 
been taken arbitrary. 

The cost function J-(x,v,t) (as well as !F{x 1 v,t)) and 
the weight factor w(x,t) describe different effects of the 
car driving, so they are not related to each other. Thus 



the terms of Eq. (|A3|) containing w, dtW, or d x w should 
be equal to zero individually. Thereby 

d x A-j t d v A = Q (A4) 

and 

A — vd v A = for w(x,t) = w a (x) , (A5a) 

d v A = for w(x, t) = w b {t) . (A5b) 

The solutions of system (|A4ll . (|A5all and system (|A4|I . 
(|A5b|) have the forms 

A(x, v, t) = C a (x) v and A(x, v, t) — Cb(t) , 

respectively, where C a (x) and Cb(t) are arbitrary func- 
tions of x and t. Thereby the desired relationship of the 
cost functions T(x,v,t), ^{x^v^t) is given for w(x 1 t) — 
w a (x) by the expression 

T{x, v, t) = xT{x, v, t) + C a (x) v (A6a) 
and for w(x,t) = Wb(t) by the expression 

F{x,v,t) = KT{x,v,t) + C b (t) . (A6b) 

where the constant a must be positive, xr > 0. 

Formulae (|A6a|) and l|A6b(l specify the admissible fam- 
ily of cost functions entering cost functional IjAlfl that 
preserve the preference relation. As seen directly, inte- 
gral IjAlfl with the integrands C a (x) and Cb(t) is inde- 
pendent of a particular form of the path x{t). So the 
transformations l|A6a|> and (|A6b|) correspond to the lin- 
ear transformation of the cost functional. 



APPENDIX B: LAGRANGIAN 
REPRESENTATION OF THE COST 
FUNCTIONAL 

The problem of finding the extremals of the functional 

/•oo 

L{h(t)}= dtw{t)T{h,v,a) (Bl) 
Jo 

defined on a set of paths {h(t)}\^ is considered. Here, 
first, the integrand T(h, v, a) is a given function of the 
headway distance h, the current velocity v = V — dh/dt, 
and the acceleration a = —d 2 h/dt 2 . Second, the weight 
factor w(t) > is a function of time t confined within 
some time interval (0,ry), i.e., w(0) ~ 1, w(t) <C 1 
for t 3> Ty, and w{t) — > sufficiently fast as t — > oo. 
So, w(0) = 1 can be adopted without lost of generality. 
Third, the trial paths {h(t)} meet the initial conditions 

h(0) = h , «(0) = vo , (B2) 

and do not diverge as time goes on, i.e., there is a number 
C so that 

\h(t)\<C. (B3) 
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The properties of the functional extremals have been 
studied already. Keeping in mind the results in Sec. 11111 
a special case is analyzed in the following. Let {r<j} be 
the characteristic time scale of the velocity relaxation to 
the stationary value V. In particular, the values of {t^} 
are estimated by the eigenvalues of the corresponding 
extremal equation linearized near the stationary point. 
The case r<j <C ry is considered here to show how the 
functional (|B1|I can be rewritten to eliminate the time 
dependent factor w(t). In this way the problem of finding 
the extremals gets the classical form met in theoretical 
mechanics. 

The stationary point {hy,V,a — 0} is determined by 
the properties of both the function T(h, v, a) and the fac- 
tor w(t) (see Exp. I|3.19|l ). As a result, at the stationary 
point (see again Exp. I|3 19|0 



d a F s t = , d h T st oc — 

TV 



(B4) 



The derivative d v T s t does not contain factors similar to 
Trf/ry. So, w(t) cannot be omitted directly. As a first 
step the function !F(h, v, a) is replaced by the difference 
AT{h,v,a) :— !F(h, v, a) — T(hy, V, 0). This does not af- 
fect the extremals. The difference A.F(/i, v,a) as a func- 
tion of t is practically confined within a time interval 
(0,Td) whose upper boundary Td > max{rd} and, so, 
Td <C Tv- Therefore, at the next step expression IjBlfl is 
rewritten as 

L{h(t)} = J^ dt^AT(h,v,a) 

-—(v-V)d v F(h v ,V,0)\ (B5) 

TV J 

which is justified to the first order in the small param- 
eter m&x{Td}/Tv by virtue of expressions (|B4|) and the 



fact that the second term in Exp. (|B5|) plays a substan- 
tial role only when the path h(t) tends to the stationary 
point. Here the symbol = means that functional l|B5|) 
possesses the same collection of the extremals as the ini- 
tial functional L{h(t)} within the adopted accuracy and 
the time dependent factor uu(t) has been approximated 
by 



w(t) 



1 



t 

TV 



(B6) 



Formula l|B6|) can be regarded actually as the definition of 
the time scale ry := [— dw(0)/dt]~ . Taking into account 
the relation (v — V) — —dh/dt and integrating the second 
term in Exp. (|B5|I by parts yields 

L{h(t)} = J^ dt^AT{h,v,a) 

- {h ~ hv) d v F(h v ,V,0)}. (B7) 

Tv > 

Omitting the constant components of the integrand 
which does not affect the extremals the required result 
follows 

L{h{t)} = J dt{F(h,v,a)- ^d v F(h v ,V,Q)}. 

(B8) 

Minimization of functional l|B8[) gives the desired ex- 
tremals and integral l|B8() does not contain a time de- 
pendent factor. The function 

F eS (h, v,a\V) := F(h, v, a) - —d v F{h v , V, 0) (B9) 

ry 

can be regarded as a Lagrangian of the car dynamics. 
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